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1  INTRODUCTION 


The  aim  of  this  work  has  been  an  investigation  into  the  distribution  of  internal  tide 
energy  over  the  Iceland-Faeroes  ridge  using  the  extensive  current  meter  data  set 
collected  by  SACLANT  and  model  predictions.  An  earlier  progress  report  (see 
appendix)  described  the  results  of  the  current  meter  analysis;  here  results  from  a 
numerical  model  will  be  compared  with  the  observations. 

Internal  tides  are  internal  waves  induced  by  the  barotropic  tide  with  energy  peaks  at 
periods  of  about  12.5  h  and  25  h.  Forcing  at  the  semi-diurnal  period  has  been 
successfully  described  by  Baines  (1982)  and  Craig  (1987)  in  terms  of  the  interaction 
of  the  flow  of  barotropic  tidal  currents  with  a  shelf  break  bottom  topography  in  the 
presence  of  oceanic  stratification.  These  models  consider  generation  in  a  two- 
dimensional  vertical  plain,  but  recent  work  by  Perkins  et  al  (1994)  and  Dale  & 
Sherwin  (1995)  has  suggested  that  in  regions  of  complex  topography  a  three- 
dimensional  approach  may  be  needed.  Thus  the  motivation  for  the  present  work  is  to 
extend  the  understanding  derived  from  two-dimensional  forcing  theory.  In  this  report 
we  look  at  the  preliminary  results  of  a  three-dimensional  primitive  equation  model  of 
internal  tidal  waves  forced  over  the  Iceland-Faeroes  ridge  by  the  output  from  a 
barotropic  tidal  model. 


2  THEORY 


The  investigations  used  the  Bryan  and  Cox  general  ocean  circulation  model,  which  is 
a  well  established  primitive  equation  model.  Only  a  brief  description  of  the  model  is 
given  here  since  it  is  extensively  discussed  elsewhere,  for  example  by  Semtner  (1986) 
and  Cox  (1984).  The  governing  equations  can  be  summarized  in  vector  notation  as: 

—+u.Vu  +  w—  +  f{kxu)  =-—Vp  +  Ay^^  +  A^V^u  (1) 

dt  dz  p  dz 


dt 


+uyT  +  w~  =  K  ^  +  K 
dz  ^  dz^  H 


(2) 


(3) 
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V7  <5h'  _ 

V-K  +—  =  0 
oz 
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where  u  is  the  horizontal  velocity  vector;  w  is  vertical  velocity;  p  is  pressure;  Ay,  Ajj, 
Ky  and  Kjj  are  the  vertical  and  horizontal  components  of  viscosity  and  diffusion;  and 
T  is  temperature.  Density,  p,  is  a  function  of  T  and  salinity,  although  for  clarity  the 
salinity  equation  equivalent  to  (2)  is  omitted  in  this  discussion.  The  differential 
operator,  V,  applies  only  to  the  horizontal  plane;  z  is  positive  upwards;  k  is  the  unit 
vertical  vector;  and  t  is  time.  Standard  simplifications  have  been  made  in  deriving 
these  equations  -  the  Boussinesq  approximation  in  (1);  the  hydrostatic  assumption  in 
(3);  incompressibility  in  (4);  and  Fickian  parametrizations  of  viscosity  and  diffusion  in 
(1)  and  (2).  For  computational  purposes  (1)  to  (4)  are  transformed  onto  spherical  co¬ 
ordinates  and  solved  on  a  finite  difference  grid. 

An  important  feature  of  the  model  is  that  is  uses  a  rigid  lid  to  suppress  the  propagation 
of  surface  gravity  waves.  Baro tropic  flows  are  represented  by  a  stream  function,  \\f, 
defined  as 


where 


{kx'V)y/=uh 
u  =  —\u  dz 


(5) 

(6) 


is  the  barotropic  velocity  vector  in  water  of  depth  h. 


(w  +  wj- 


dr 

The  internal  tide  is  forced  in  the  present  configuration  by  replacing  w  —  in  (2)  with 

dz 

—  where  Wq  is  the  oscillating  vertical  velocity  due  to  the  barotropic  tide. 
az 

This  is  equivalent  to  the  linear  forcing  technique  used  by  Baines  (1982),  although  it  is 

incomplete  because  it  omits  terms  such  as  (u-i-Uo)-VT  where  Uq  is  the  horizontal 

cfr 

current  vector  of  the  barotropic  tide.  However,  except  near  a  front,  we  expect 
to  be  much  greater  than  Uq-VT. 


3  DERIVATION  OF  BAROTROPIC  FORCING 

The  barotropic  forcing  term,  Wq,  has  been  derived  from  a  barotropic  model  of  the  tides 
of  the  NE  Atlantic  due  to  Roger  Hather  at  the  Proudman  Oceanographic  Lab.  The 
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model  results  specify  the  amplitude  and  phase  of  vertical  elevation  (which  is  ignored) 
and  horizontal  velocity  for  different  tidal  constituents.  In  the  presence  of  variable 
topography  the  continuity  equation  of  the  barotropic  tide  can  be  written  as 


where 


^+v.(e/!-')=o 


Q-hiu  cosicn-G  )i+v  cosiot-G^^  )j) 
^  V  (7  ^  MCF''  a  ''  Va'^ ' 


(7) 

(8) 


is  the  volume  flow.  Here  a  is  the  frequency  of  a  specific  tidal  constituent;  u^j  and 
are  the  amplitude  of  the  barotropic  tidal  currents  in  the  i  and  j  directions;  and  Gy  and 
Gy  are  the  corresponding  phase  lags.  On  the  assumption  that  V-Q  is  small  the 
expression 


.^  =  -l2-vr'  (9) 

dz 

has  been  used  to  derive  Wq.  In  the  Cox  model  the  expression 

yv^{z)  =  -Q-Vhr^z  (10) 

has  been  used  to  give  the  vertical  velocity  at  depth  z  with  the  rigid  lid  assumption. 
The  computation  of  Wq,  is  complicated  by  the  fact  that  the  Flather  model  uses  an 
Arakawa  'C  grid,  whilst  the  Cox  code  uses  an  Arakawa  'B'  grid. 


4  THE  GRIFFITHS  MODEL  OF  THE  ICELAND-FAEROES  RIDGE 

The  Griffiths  (1995)  modification  of  the  Cox  code  rigid  lid  General  Ocean  Circulation 
Model  has  been  used  to  establish  the  background  density  field  for  the  Iceland-Faeroes 
ridge  region.  The  model  parameters  are  given  in  Table  1.  The  topography  was  taken 
from  the  DBDB5  data  set  and  the  model  initialized  with  a  temperature  and  salinity 
field  from  the  Levitus  data  set.  The  model  computes  baroclinic  fields  of  density  and 
momentum  and  a  separate  barotropic  momentum  (or  stream  function)  field.  Passive 
boundary  conditions  are  applied  to  the  baroclinic  fields  and  to  the  north  and  south 
stream  function  boundaries.  However,  values  of  the  stream  function  on  the  east  and 
west  boundaries  are  forced  to  equal  values  predicted  by  a  larger  scale  model  due  to 
Stevens  (1990).  The  model  was  run  for  typically  five  days  during  which  time  inertial 
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oscillations  (period  about  13.7  h),  which  are  caused  by  the  initialization  process,  were 
gradually  damped  out. 


5  TIDAL  FORCING  AND  ANALYSIS 

Throughout  a  model  run  the  internal  tide  was  forced  using  (10)  over  the  model 
domain.  For  the  last  wave  period  of  a  run  (duration  T  =  27r/G)  an  analysis  of  the  tidal 
current  in  the  i  direction  at  each  grid  point  was  conducted  using 

U{x,y,  z)  =  -r{^u{x,  y,z,  t)  exp(iOnA  J  (11) 

72=1 

where  U  is  complex;  u  is  the  speed  of  the  current  at  time  t  in  the  i  direction  and 
(NAj  =  T).  A  similar  analysis  was  conducted  to  derive  V  with  v  (the  speed  in  the  j 
direction)  instead  of  u  in  (1 1).  If 

U{x,y)  =  ^\  U{x,y,z)dz  (12) 

is  the  i  component  of  the  (complex)  depth  mean  current,  and  V  is  the  equivalent  j 
component,  then  the  baroclinic  components  of  the  tidal  current  are 

Uf,{x,y,z)^U{x,y,z)-U{x,y)  (13a) 

and 

Vj(x,y,z)  =  y(x,y,z)-y(x,y)  (13b) 

A  baroclinic  tidal  energy  parameter  was  thus  computed  as 

*)dz  (14) 

where  u^,*,  Vj,*  are  the  complex  conjugates  of  u^  and  Vj^. 


6  RESULTS 


The  model  was  separately  forced  with  two  tidal  constituents,  M2  (period  12.42  h)  and 
01  (period  25.82  h)  using  Levitus  data  to  prescribe  a  climatic  spring  time  density 
field.  The  distributions  of  r  for  each  constituent  within  a  box  bounded  by  67°N  15°W, 
62°N  3°W  (within  the  Griffiths  model  domain)  are  compared  in  Figs  1  and  2  with 
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similar  plots  derived  from  current  meter  observations  (described  in  Progress  Report 
No  2,  see  Appendix). 


The  predictions  with  M2  are  about  a  factor  of  4  smaller  than  those  observed,  but 
within  the  comparison  area  encouragement  can  be  drawn  from  the  fact  that  the  two 
energy  distribution  patterns  seem  comparable.  For  example,  large  values  are  both 
predicted  and  observed  along  the  SE  edge  of  the  Iceland  shelf,  whilst  they  are  much 
smaller  along  NE  shelf  edge.  Overall  we  notice  that  the  model  predicts  relatively  high 
semi-diurnal  internal  tide  levels  along  the  crest  of  the  ridge  which  is  reasonable  since 
strong  barotropic  tidal  currents  cross  it. 


These  patterns  are  not  reproduced  with  01 .  Here  the  magnitude  of  the  predictions  and 
observations  are  very  similar  on  the  SE  edge  of  the  Iceland  shelf,  but  the  predictions 
are  very  much  smaller  than  the  observations  on  the  north  side  of  the  ridge  and  on  the 
NE  Icelandic  shelf  edge. 


7  CONCLUSIONS 


This  work  has  demonstrated  that  it  is  worthwhile  trying  to  map  internal  tide  energy 
levels  in  regions  of  complex  topography.  The  method  of  coupling  the  baroclinic 
mode  of  a  rigid  lid  general  ocean  circulation  model  to  forcing  from  a  barotropic  tidal 
model  is  stable  and  permits  runs  to  take  place  in  a  realistic  time  (the  results  shown 
here  took  about  2  h  on  a  Sun  Sparc  10).  However,  comparison  with  observations 
shows  that  further  work  is  required  before  the  results  can  be  treated  with  confidence. 
For  example  it  will  be  necessary  to  look  at  the  sensitivity  of  the  predictions  to  the 
prescribed  density  field.  It  is  likely  that  the  stratification  used  to  the  north  of  the 
Iceland-Faeroes  ridge  is  much  weaker  than  actually  exists  and  may  explain,  in 
particular,  why  the  predicted  01  magnitudes  are  too  small  in  that  region.  It  may  also 
be  worthwhile  to  examine  the  sensitivity  to  vertical  resolution  in  the  surface  layers. 


It  should  be  evident  that  this  contract  has  stimulated  some  innovative  work  which  is 
not  yet  complete.  This  work  will  be  continued  with  funding  from  another  agency  and 
will  eventually  be  reported  in  the  scientific  press  with  due  acknowledgement  to  the 
support  of  NRL. 
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Table  1 


The  Model  Parameters 


NW  comer 

67°N  20°W 

SE  comer 

59°N  0°W 

Number  of  grid  points  (xxyxz) 

121x97x11 

Ax  (longitudinal  grid  spacing) 

5’ 

Ay  (latitudinal  grid  spacing) 

10’ 

A2  (vertical  grid  spacing,  from  surface) 

50,  50,  50,  100,  100,  100,  100, 
100,  300,  300, 1500  m 

Af  (time  interval) 

900  s 

Ah 

5x10®  cm2  §-1 

Av 

1  cm2  s-i 

Kh 

1x10®  cm2  s-i 

Kv _ 

1  cm2  s-i 
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Estimates  of  r  for  the  internal  tidal  constituent,  M2,  in  cms'* 
a)  from  observations;  b)  from  model  predictions. 


LONGITUDE 


Longitude  (deg  W) 


Estimates  of  r  for  the  internal  tidal  constituent,  01,  in  cms'^ 
a)  from  observations;  b)  from  model  predictions. 


APPENDIX 


To  Model  Tidal  Structure  in  the  Shallow  Water  of  the  Iceland-Faeroes 

Ridge  Region 

Progress  Rpt  No  2 

Project  UCES  105 

ONR  Grant  NOOO 14-93- 1-60 11 

Dr  TJ.  Sherwin,  UCES,  Marine  Science  Labs,  Menai  Bridge,  Gwynedd,  UK 


1.  Introduction 

The  main  objective  of  this  work  is  to  investigate  the  distribution  of  internal  tide  energy  over 
the  Iceland-Faeroes  ridge,  using  (i)  the  extensive  existing  current  meter  data  set  available  and 
(ii)  predictions  from  the  UCES  two-dimensional  linear  slice  model  of  internal  tide 
generation.  To  date  we  are  close  to  finishing  the  first  part  of  the  work  and  this  report  gives  a 
summary  of  the  observed  distribution  of  two  semi-diurnal  constituents  (M2  and  S2)  and  two 
diurnal  constituents  (K1  and  01).  We  are  about  to  start  the  second  stage  of  the  work  which 
involves  explaining  the  distribution  using  the  linear  model.  The  final  phase  of  the  work  will 
be  the  preparation  of  a  report,  the  main  part  of  which  will  be  published  in  a  scientific  journal. 


2.  Data  Acquisition 

The  observations  described  below  all  come  from  the  SACLANT  archive.  We  have  trawled 
the  relevant  data  banks  (principally  the  British  Oceanographic  Data  Centre)  for  archived 
current  meter  data  of  the  region.  In  the  event  only  a  small  number  of  useful  current  meter 
sites  have  been  located  -  in  the  archives  of  the  Bundesant  fiir  Seeschiffahrt  und 
Hydrographie,  Hamburg.  These  observations  were  made  in  1973  using  Richardson  current 
meters  and  are  of  variable  quality.  So  far  they  have  not  been  included  in  the  analysis 
although  they  cover  the  south  eastern  half  of  the  ridge  for  which  there  is  no  other  data. 


3.  Data  Evaluation 

The  SACLANT  current  meter  data  archive  have  been  analysed  by  our  partners,  under  Dr  H 
Perkins,  and  his  results  are  discussed  below.  The  method  of  analysis  was  as  follows:  depth 
integrated  (barotropic)  currents  were  computed  at  each  site  and  compared  with  the 
predictions  from  the  tidal  model  due  to  R.A.  Flather.  From  this  we  derived  sufficient 
confidence  in  the  Flather  model  to  use  it  as  a  reliable  estimate  of  barotropic  tidal  currents 
over  the  whole  domain  (Fig.  1).  This  was  necessary  because  on  some  rigs  there  existed  too 
few  current  meters  to  justify  the  computation  of  a  depth  mean.  The  Flather  tidal  currents 
were  then  subtracted  from  the  individual  current  meter  time  series  and  a  tidal  analysis  was 
conducted  on  these  residual  time  series  to  produce  individual  estimates  of  the  internal  tide  for 
each  of  the  constituents,  M2,  S2, 01  and  Kl. 

From  these  constituents  estimates  of  the  internal  tide  activity,  r,  have  been  calculated  as 
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2N 


+  V 
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where  N  is  the  number  of  current  meters  on  a  rig,  and  u  and  v  are  the  amplitude  of  the 
internal  tidal  current  in  the  east  and  north  directions  respectively.  Values  of  r  are  shown  in 
Figs  2  to  5. 

There  is  a  clear  difference  in  the  spatial  distribution  of  activity  between  the  super-inertial, 
semi-diurnal  constituents  (M2  and  S2)  and  the  sub-inertial,  diurnal  ones  (K1  and  Ol).  The 
peaks  in  activity  of  M2  and  S2  are  found  in  the  region  where  the  ridge  meets  the  Icelandic 
shelf  and  to  some  extent  along  both  flanks  of  the  ridge  on  its  shallower  northern  end  (Figs  2 
and  3).  By  contrast  K1  has  relatively  high  values  around  the  Icelandic  shelf,  with  some 
outlying  regions  of  higher  activity  on  both  flanks  of  the  ridge  (Fig.  3);  and  the  Ol  internal 
tide  is  mainly  located  on  the  eastern  side  of  Iceland  and  the  northern  flank  of  the  ridge 
(Fig-  4). 

There  are  probably  two  explanations  for  this  variability  which  will  be  investigated  in  the 
next  part  of  this  work.  One  reason  is  that  the  spatial  distribution  of  the  barotropic  tide 
(which  is  the  forcing  agent)  appears  to  vary  in  a  way  which  is  compatible  with  the  internal 
tide  observations  (see  Fig.  1).  It  is  particularly  noticeable  that  the  internal  tide  distributions 
of  K1  and  Ol  match  their  barotropic  counterparts  quite  well.  The  other  explanation  concerns 
the  fact  that  on  a  rotating  Earth  super-inertial  waves  can  propagate  into  the  ocean  away  from 
their  source  whereas  sub-inertial  ones  remain  trapped  to  the  topography.  Perkins  et  al. 
(J.P.O.,  1994)  have  shown  that  the  internal  tide  generated  at  the  edge  of  the  Iceland  shelf  has 
the  potential  to  propagate  a  significant  way  across  the  ridge. 


4.  Future  Work 

It  is  planned  to  compliment  the  analysis  of  the  observation  with  a  series  of  two-dimensional 
internal  tide  model  runs  along  sections  that  will  make  it  possible  to  map  the  distribution  over 
most  of  the  region  shown  in  the  figures.  If  these  mapped  predictions  agree  with  the 
distribution  seen  in  the  observations  then  we  will  have  confidence  in  the  estimates  thus 
produced  and  this  may  provide  the  stimulus  to  conduct  a  similar  analysis  in  other  shelf  edge 
and  shelf  regions. 
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Figure  1.  The  amplitude  of  the  eastward  component  of  current  velocity  for  the  four  major 
tidal  constituents  over  the  Iceland-Faeroes  ridge  (Kl,  01,  M2  and  S2)  from  the  model  due  to 
R.A.  Flather.  Units  are  in  cms  s-‘,  but  note  that  the  contour  intervals  vary  between  the 
figures.  The  northward  currents  show  a  similar  distribution. 
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Figure  2.  Estimates  of  r  for  the  internal  tidal  constituent,  M2.  Units  are  in  cms  S'^. 
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Estimates  of  r  for  the  internal  tidal  constituent,  S2.  Units  are  in  cms  S‘h 
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Figure  4.  Estimates  of  r  for  the  internal  tidal  constituent,  Kl.  Units  are  in  cms  s  h 
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Figure  5.  Estimates  of  r  for  the  internal  tidal  constituent,  01.  Units  are  in  cms  s  h 
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